NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 1942 


THE EFFECT OF WALL FRICTION ON THE STRENGTH OF SHOCK 
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SUMMARY 


A theory is presented for the attenuation due to viscous boundary 
layer of shock waves traveling in closed passages. Because of the 
similar nature of the problem, the theory is also developed for the 
attenuation of a hydraulic Jump moving in a channel . The results of 
experiments that bear out the theory are presented for both a shock 
tube and a hydraulic channel. 

INTRODUCTION 


Recent interest in intermittent— flow devices, as, for example, 
shock tubes, has caused questions to be raised concerning the 
viscous effects. 


Although a large amount of literature exists on the theory of 
shock tubes in which perfect fluids are used (see, for example, refer- 
ences 1 and 2), no discussion is known to be available of the effects 
of viscosity in such devices. The present analysis is an attempt to 
develop a quantitative method for determining the attenuation of a 
shock wave moving in a tube or passage resulting from the action of 
viscosity at the walls. 

Because of the similar nature of the attenuation of a hydraulic 
Jump traveling in a channel, the method is also developed for this 
phenomenon . 


SYMBOLS 


a velocity of small disturbances in fluid 


A 


area of tube 
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g acceleration due to gravity 

h depth of liquid 

H depth ratio (h 2 /h-[) 

z(h) = (H - i)\/h(h + 1) 

m mass flow per unit time 
M Mach number 

n(F) = \J6? + 1 

p pressure 

P pressure ratio (P 2 /P 1 ) 

r radial coordinate from center of tube 

E radius of cylindrical tube 
t time 

T absolute temperature 
u velocity 

w width of channel 

x horizontal coordinate 
y vertical coordinate 

7 ratio of specific heats 
H absolute viscosity 

v kinematic' viscosity (n/p) 

p density 

Subscripts: 

0 conditions in fluid on high-pressure side of diaphragm or barrier 

1 conditions in fluid on low-pressure side of diaphragm or barrier 
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2 conditions in fluid "behind- shock or hydraulic Jump 

3 conditions in fluid after wave reflection 

ex pertaining to expansion wave 

i pertaining to initial flow conditions 

s pertaining to shock or hydraulic Jump 

BOUNDART LATER IN A TUBE IN WHICH THE ENTIRE 
FLUID IS SUDDENIT ACCELERATED 


Since the "boundary layers that cause the attenuation considered 
are assumed to depend principally on time, the differential equation 
of viscous motion is solved for the case of an infinitely long tube 
of radius R in which the entire fluid at time t = 0 is suddenly 
set in motion at velocity u^. In this case, the velocities caused 
by viscous action depend only upon the distance from the center of 
the tube and on the time that the fluid has been in motion, rather 
than upon position along the length of the tube; therefore, the 
differential equation of viscous motion in cylindrical coordinates 
becomes 


du 

^t 



( 1 ) 


The boundary conditions to be imposed are: 


at t = 0 


and at r = R 


u = Uj^ 


(0 < r < R) 


> ( 2 ) 


u = 0 

fn addition to these boundary conditions, the function 
be continuous at r = 0. 


(t > 0) 


J 

u = f(r,t) must 
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The solution of this equation subject to these bpundary conditions 
is well-known in the theory of heat conduction (see reference 3) and is 


u 

u i 


T - a a 2 ' r tM i 

2 f J °v~> 

n=l anJi(an) 


( 3 ) 


where the values of oq are the roots of the equation Jo(x) = 0 
and Jq and J]_ are Bessel functions of the first kind of order 0 

and 1, respectively. 

The velocity profiles in the tube for various values of the time 
parameter vt/R2 are plotted in figure 1. This solution for sma ll 
values of the time parameter reduces to the two-dimensional solution 


where 


u 

u i 



erf (x) 


_2_ 

ft 



dx 


and y is the distance from the surface at which the velocity is 0 
and in this case is equal to R — r. This solution is useful not only 
for evaluating the boundary layer in tubes of other than circular cross 
section but also for evaluating the series of equation ('3) for small 
values of the time parameter since, for these values of the parameter, 
many terms of the series are necessary.. 

Figure 1 indicates that, as time passes, the mass flow in the tube 
decreases from its maximum and initial value at t = 0 of m* = npuiR^, 


The ratio of the mass flow at any time to the initial mass flow 

Jo 2 ’"™ 6 * 

• m 


00 

= V — 2 0 

jtpUiR^ n=1 an 


00 


is plotted in figure 2 as a function of the time parameter vt/lf 
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For small values of the time parameter equation (4) reduces to the two- 
dimensional solution 



where A is the area of the tube, so that for the circ ular tube 



This result is also plotted in figure 2 for comparison with the results 
of equation (4). These solutions are used to discuss the attenuation 
of a shock wave traveling in a tube. 


ATTENUATION OF A SHOCK WAVE IN A TUBE 


A shock tube is a closed tube that is divided into two compartments 
by means of a breakable diaphragm. If the two compartments are pumped 
to different pressures p 0 and pj_ where p 0 > Pp and the diaphragm 

is broken, an expansion wave travels into the 'compartment where the air 
is at rest at p 0 j the air is thus accelerated to velocity Ug and the 
pressure decreased from p 0 to p 2 . A shock wave travels into the 
compartment where the air is at rest at pressure p^_, accelerates 
it to the same velocity U2, and raises its pressure to the 
same pressure P 2 - This process is illustrated in figure 3. The result 
of bursting the diaphragm in the absence of viscosity is to set the 
entire mass of air between the expansion wave and the shock wave in • 
motion at velocity u 2 . Thus, the mass flow per unit t ime set in 
motion across a section of the tube between the diaphragm and the 
shock wave is 


m — np 2 u 2 p2 
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This mass flow is related to the pressure ratio across th9 shock 
Po 

wave — = p and the conditions ahead of the wave by the following 
p A 

equation (for 7 = 1.4) 


m 


np-ja-jR 2 ff 


= -^vEpTI ~ 1 = 5- n(p) 


P + 6 




(5) 


which is' derived in the appendix; thus. 


n(P) « \/6P + 1 


P - 1 
P + 6 


( 6 ) 


In order to evaluate the effect of viscosity upon the flow in the 
tube, the following assumption is made: The thickness of the boundary 

layer at any point behind the shock depends only upon the length of time 
the air has been in motion with respect to the wall of the tube, or upon 
the length of time since the shock wave passed those particles. Such an 
assumption neglects variations of velocity in the stream direction and 
can be shown to give an answer to the steady flat— plate boundary-layer 
flow that Is a good approximation to the Blasius solution. Since, in 
the present problem, the shock wave travels faster than the fluid in 
motion behind it, this assumption should be an even better approximation 
to the real solution. Thus, the boundary layer between the shock wave 
and the particles that were originally at the diaphragm is essentially 
as shown in figure 4.' 

When this assumption is made, the solution obtained In the 
preceding section is applicable; If attention is confined to the 
position in the tube at which the particles, originally at the diaphragm 
are found, the mass flow in this section must drop off with time according 
to equation (4), so that the mass flow at any time t is 



n=l 



-an S vt/E 2 


m = m± 4 
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If the value of the initial mass flow in terms of the initial 
pressure ratio (equation ( 5 )) is substituted in this expression, the 
mass flow at any time is given by 


m 


np-ja-^ 




P< - 


-1 1. T -L 


( 7 ) 


This formula indicates that the mass flow in the tube decreases with 
time according to the function that is plotted in figure 2. 

Furthermore, the shock wave itself sets in motion a mass flow 
that is given by (see equation ( 5 )) 


m 

jrp^jR 2 


= x/SFTI 
V 7 


p - 1 
P + 6 


If this mass flow is greater than that crossing the tube at the 
position of the diaphragm particles, an expansion wave would be formed 
that would catch up with the shock wave, reduce its . strength, and thus 
reduce the mass flow it was initiating; but if the mass flow caused by 
the shock was less than that across the diaphragm— particle section, a 
compression wave would develop that would move after the wave, catch it, 
increase its strength, and thus increase the mass flow it was initiating. 
Hence, it may be assumed for the purposes of estimating the attenuation 
that the mass flow initiated by the wave is just equal to that given by 
equation ( 7 ). If both the expansion wave and the compression wave were 
reduced in this manner, a pressure gradient would be produced between 
these two waves in order to preserve the pressure ratio between the two 
undisturbed regions. Part of this pressure gradient would be used to 
accelerate somewhat the flow between the two waves, and thus the decrease 
in mass flow at the diaphragm— particle section would be somewhat less 
than that given by equation ( 7 ). However, for the purpose of determining 
the approximate magnitude of the attenuation, the effect of this pressure 
gradient is neglected. If the mass flows given by equations (5) and (7) 
are set equal to each other, the relationship between the strength of 
the wave P and its Initial strength Pj is found to be 


\J6t + 1 


P - l 
P + 6 


= v^Pi + 1 


Pi 


1 -an-vt/BP 


( 8 ) 


n=l an 
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The function 
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n(P) = \jS P + .1 


P - 1 
P + 6 


is plotted in figure 5 as a function of the strength of the wave P. 
Substituting n(P) from < equation (6) in equation (8) gives 

n(P) = 4 _l_ e -% 2vt / E2 (9 

n(Pi) n=l cm 2 

which is the same function as is plotted in figure 2. 

If the strength of the shock at any distance from the diaphragm is 
to he determined, the values of t need only he found in terms of the 
distance x and the velocity of the shock u s and the kinematic 
viscosity behind the shock evaluated in terms of the shock strength P 
and the conditions ahead of the shock. If the attenuation is small. 



Then, since u„ is a function of the strength of the shock and the 
s i 

condition ahead of the shock, V 2 t/R 2 becomes (see derivation in the 
appendix) 


Vgt = /F(Pi + 6)^ Vl x 
R 2 (6P* + l) 2 * 5 a-^ 2 


( 10 ) 


where linear variation of viscosity with temperature is assumed. 

If the initial strength of the wave and the conditions into which 
it is traveling are known, v 2 t/R 2 can be evaluated for any distance 
from the diaphra,gm. After this value is determined, the ratio n(p)/n(P i ) 
may be found from figure 2. Then, from figure 5, the value of the 
strength of the wave P . at x may be evaluated. 

Before proceeding with an account of an experimental test of this 
theory, the theory for the attenuation of a hydraulic jump in a channel 
is presented. 
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BOUNDARY LAYER IN A FLUID OF FINITE DEPTH 


An analysis parallel to that of the boundary layer in a cylinder 
may be carried out for the boundary layer in a fluid of finite depth on 
a flat plate. This analysis can be applied to find the attenuation of 
a hydraulic jump traveling in a channel. 

The differential equation (corresponding to equation (l)) is 

du _ v cj£u 

^t cJy2 


The boundary conditions are: 
at t = 0 

at y = 0 


(0 < y h) 


at y = h 


u = 0 (t > 0) 

- 0 (t > 0) 


The last condition expresses the absence of tangential force on the 
free surface. 


The solution of this boundary— value problem is also known from 
the theory of heat transfer. (See reference 3- ) The solution is 


u 

“I 



sin 



* h 2 


The ratio of the mass flow past, a section of width w at any time t 
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to its value at t = 0 of 


m^ = wphu^ 


is 



(an - 1)2 


-#( 2 n-l ) 2 „2 


Vt 

h 2 


( 11 ) 


Tliis result is plotted in figure 6 and is used to evaluate the attenu- 
ation of a hydraulic Jump traveling in a channel. 


ATTENUATION OF A HTDEAULIC JUMP IN A CHANNEL 


If a channel is divided into two parts by a barrier and the two 
parts are filled with a liquid to different levels h Q and hq, 
where ho > hq, and the barrier is suddenly removed, a situation 
analogous to that of a shock tube exists. Although the equations of 
a shock and a hydraulic Jump are completely different (not Just the 
same equation using different ratios of specific heats, as in the 
case of continuous flows such as the water analogy, see reference h) 
the phenomena are fundamentally similar ; in particular, the liquid 
originally at a depth hj_ is suddenly accelerated to a velocity u^ 
and the depth is increased to hp when the Jump passes by it. 

t 

The mass flow per unit time set in motion immediately behind a 
hydraulic Jump in a channel of width w is 


m = pwu/^hg 

This mass flow is related to the ratio of the depths on the two sides 
ho 

of the Jump — = H and to the conditions ahead of the wave by the 
h l 

equation (derived in the appendix) : 


m 

pwa-jh^ 


= ~ (H — 1)\/H(H + 1) = -L Z(H) 

V2 V2 


( 12 ) 
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where a.± = \Jgh-j_ is the speed of an infinitesimal ware. The 

parameter 1 (h) is plotted as a function of the depth ratio H 
in figure T« 

By an argument parallel > to that for a shock tube, if the channel 
is assumed to he wide in comparison with the depth of the liquid so 
that edge effects may he neglected, the strength of the wave after it 
has gone a distance x can he determined hy the relation 


i(h) 8_y 

l(Hi) , n 2 n^l 


1 

(2n - l) 2 



(2n— l)2-^£- 

h.2 2 


(13) 


If vt/hg 2 is expressed in terms of the Jump strength and condi- 
tions ahead of the wave, and if the assumption is made that t = 

U Si 

this parameter hecomes, as shown in the appendix: 


vt _ _1_ / 2 yi 

h 2 2 " H 2 VH(H - 1) ^2 


(l4) 


With the aid of equation (l4) and figures 6 and 7, the strength of a 
hydraulic Jump at a position in the channel can he computed. 


TESTS OF THE ATTENUATION OF A SHOCK WAVE DT A TUBE 


In order to obtain an experimental check of the theory Just 
presented, a special shock tube was constructed. The tube was 0.4 inch 

in diameter, and the high-pressure compartment was 4^ feet from the 

diaphragm to the end of the tube. Two low-pressure chambers were 
interchangeable on the other side of the diaphragm — the first 1 foot 
long, the other 11. 3 feet long from the diaphragm to the pickup 
crystal. The experimental apparatus with the short low-pressure 
chamber is shown in figures 8 and 9- 

Tests were made by charging the high-pressure chamber with nitrogen 
to the desired pressure. The diaphragm was then broken by puncturing it 
with a solenoid-driven plunger, and a shock wave passed into the 
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low— pressure chamber, which was at atmospheric pressure, reached the 
end of the tube, and was reflected back toward the diaphragm. The 
* voltage output of a piezoelectric crystal mounted flush in the end of 
the tube was recorded. This voltage, which was a measure of the ratio 
of the initial pressure in the tube and the pressure after the shock 
wave had been reflected, was recorded by means of an oscillograph and 
rotating drum camera the same as that described in reference 2. One 
such record, which is typical, is presented herein as figure 10. 

By this method the voltage output of the crystal was obtained for 
a range of initial ratios of pressure between the high-pressure and 
low-pressure sides of the diaphragm with both the 1— foot and the 
11.3— foot pressure chambers. Since the reflected— pressure ratio 
measured (p^/pp) is related to the strength of the wave striking the 

end of the tube by the formula (see the appendix) , 

^3 _ P(8P - l) 
p l P + 6 


the records made are a measure of the relative attenuation of the wave 
at stations 1 foot and 11.3 feet from the diaphragm said are seen to be 
independent of the crystal characteristics. 

The theoretical reflected— pressure ratio p^/Pi &s a function of 
initial— pressure ratio p^/pp, with and without the effect of viscosity 

considered, is plotted in figure 11 for a station 1 foot from the 
diaphragm, along with the voltage output of the crystal measured at 
the same station. The derivation of the relation between Po/Pi 

and P is given in the appendix . 

If the theory is assumed to be correct, these data may be used to 
determine the voltage output of the crystal for a given pressure rise. 
This information is given in figure 12. The reflected— pressure ratios 
corresponding to the crystal voltage outputs measured with the 11. 3— foot 
tube can be determined from this figure and are plotted in figure 13, 
along with the theoretical pressure ratio at this station. 

Although the theory seems always to predict too little attenuation, 
the order of magnitude (12 to 15 percent) agrees very well with that 
measured (12. 5 to 16.7 percent). However, the theory takes into account 
only one factor contributing to the attenuation, and other factors (such 
as the impossibility of a perfectly efficient diaphragm) may contribute 
to a reduction of the strength of the initial wave formed which might 
well be the order of the difference between the theory and experiment. 
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TESTS OF THE ATTENUATION OF A HYDRAULIC JUMP IN A CHANNEL 


For experimentally determining the attenuation of hydraulic Jumps 
an apparatus was used that consisted of a standard 8-inch steel channel 
20 feet long with plates welded to its ends. A transverse cut was made 
in the middle, into which a removable metal barrier was inserted. Liquid 
was poured in to different depths on the two sides of the barrier. The 
barrier was suddenly removed, and the depth Just behind the wave was 
measured at two points — one close to the barrier, and one farther down- 
stream. The depths were measured in sixty— fourths of an inch by means of 
thin steel scales. By this means the height of the moving liquid could 
be obtained within 3 percent. 

The first tests were made with water. The attenuation of the waves 
was found to be about twice that predicted by the theory. However, 
surface tension, unaccounted for in the theory, was believed to contrib- 
ute a large share of the attenuation. In order to avoid this effect, 

SAE 10 oil, which has a surface tension about one— third that of water 
and a kinematic viscosity about seventy times that of water, was used 
bo that the effect of surface tension would be negligible compared to 
that of the viscous boundary layer. Measurements were made of the 
attenuation of hydraulic Jumps of various strengths moving into liquid 
depths of from 30/64 inch to 60/64 inch. The height measurements for 
this purpose were made at a position 2 inches downstream and at a 
position 4 feet downstream from the barrier. 

Some of the results are given in figure 14. Since the derivation 
of the relation between Hj_ and ho/hq is straightforward (see 
appendix) the experimental values of the hydraulic Jump strength H 
measured close to the barrier position plotted against ho/hq when 

compared with the theoretical curve indicate the accuracy of the 
measurements . On the other hand, the plots of H measured 4 feet 
downstream against ho/hq indicate the reliability of the theory 
advanced herein. In the case of the initial strengths an error of 
measurement seems to exist which gives wave strengths consistently 
higher than the theoretical values. If the results are corrected 
for this consistent error, excellent agreement between theory an d 
experiment is obtained. 


CONCLUSIONS 


A theory has been presented for the attenuation due to viscous 
boundary layer of shock waves traveling in closed passages. Because 
of the similar nature of the problem, the theory was also developed 
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for the attenuation of a hydraulic jump moving in a channel. The theory 
in its two forms was then checked experimentally. 

Experimental results obtained by two greatly different methods 
indicate that the the’ory is sufficiently accurate to allow analysis of 
the effects of viscosity in shock tubes and intermittent— flow apparatus. 

The assumption that the boundary layer in the intermittent— flow 
processes described depends only on time seems to be Justified. 

Experiments made in the shock tube show that, in the relatively 
small tube used, the pressure ratio across the shock wave dropped only 
some 12 to 15 percent in a distance of 11 feet. Thus, since the parameter 
controlling the attenuation is increased in proportion to length and to 
the inverse square of the radius, in the passages of a normal intermittent- 
flow apparatus the effect of viscosity at the walls causes little change 
in the strength of the waves moving in the passages. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., July 6, 1949 
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APPENDIX 

DERIVATION OF FORMULAS 
Derivation of Formula for m/np^ajR 2 


The equations characterizing a moving plane shock may he found by 
changing the frame of reference so that the shock appears to be 
stationary. The velocities in this frame of reference are denoted by 
primes; the other quantities remain unchanged in the two frames of 
reference. The equations of a stationary normal shock are well-known . 
(See reference 5*) Thus, if the ratio of specific heats 7 is taken 
as 1.4, the pressure ratio across the shock is given by 



( 15 ) 


and the density ratio by 


where 


P 2 _ 6 M ]_ 2 

P 1 "V + 5 


(16) 



(IT) 


Further, the continuity equation is 


Up’pp = U^Pg 


(18) 


In order to change the frame of reference back to the case of a shock 
moving into a still gas, the relations 



j 
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are used. By substituting values from equations (15) to (19) in the 
expression for mass flow 


m = itPpUpR 2 

the following relation, equation (5) in the text;, is obtained: 


2 = -p= 

Itp i a i^' ft p + 6 


Derivation of Formula for Vpt /H 2 
If the gas law is written 


P = !a = ^2 T2 

pi pi t-l 


and if viscosity is assumed to vary linearly with temperature 


P2 _ Tg 
“1 = T i 


the kinematic viscosity can he written 


M-p Ui / Pi \ 

Vp = — = — p( — ) = V.,P| 

Pp P^_ VPp/ 


V P 2y 


( 20 ) 


Writing t = — and using equations (15) to ( 20 ) result in 
u s 

v 0 t 


v 2 t _ V7(P + 6) 2 P v l x 

R 2 (6P + l) 2 '^ a^R 2 


which is equation (10) in the text. 
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Derivation of Formula for Pq/Pi 

Reference 6 shows that, on the two sides of an expansion wave, the 
relation 



( 21 ) 


holds. In this case, Uq = 0 and Sq = ap since the original tempera- 
ture on the two sides of the diaphragm is the same. If the speed of 

sound is taken as a = . /^£, the isentropic relation as — = 

VP Pi 

7 as lA in equation (2l), and by use of equations (15) to (19), the 
following expression is found for Pq/Pi : 




P -1 1 7 

ff(6p + i)_ 


Derivation of Formula for p^/pp 

In order to determine the pressure in the end of the tube after the 
wave is reflected, a new frame of reference is taken so that the 
reflected wave appears to be stationary. This change is accomplished by 
setting 


u l’ 

u 2* 






where u a _ is the velocity of the reflected shock. 


( 22 ) 


Equations (15) to ( 19 ) can be solved for U2 in terms of P, but 
equations (15) to (l8) also apply to the reflected shock (after suitable 
changes in the subscripts) and, together with equation (22), can be 
solved for ug in terms of P^/Pg* Equating these two expressions 
for u 2 and simplifying result in the expression 


£3 = 8p - 1 

P 2 P + 6 


£3 _ p(8p - 1) 

Pi' - P +6 


Hence 
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Equations for Hydraulic Jump 

If a frame of reference is taken so that the jump appears to he 
stationary, the equation of continuity can he written as (see 
reference 4) 

u l* h i ; = u 2 ,h 2 (23) 

and the momentum equation as 

+ U^* 2 !^ = ^ — + Ug^l^J (24) 

By substituting values from equations (19), (23), and (24) in the 
expression for mass flow 


m = pwugh^ 

the following relation, equation (12) in the text, is obtained: 


m 

pwajh-L 


- l)i/H(H + 1) 

\[2 


Also, the following expression for vt/h 2 2 ,. equation (l4) in the text, 
is found: 


Vt _ _1_ / 2 v x 

h 2 2 H 2 YH(H + l) a^ 2 


Since the expansion wave is a continuous phenomenon, equation (21 ) 
may he used in the hydraulic case in which 7 = 2 .' (See reference 4 .) . 
Thus, if the speed of a smn.1 1 wave is a = Vgh and if Uq = 0 in 
equation (2l), and by use of equations (19), (23), and (24), the following 
expression for hQ/h^ is obtained: 

^ = [/ff ^ (H - 1) 
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Figure 4.- Boundary layer on a wall behind a shock wave traveling into air at rest. 
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Figure 5.- Plot of the function n(P) from equation (8). 
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Figure 6.- Ratio of mass flow to initial mass flow in a shock wave given by 
equation (11) as a function of the time parameter . 


1(H) = (H - 1)>/H(H + 1) 
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Figure 7.- Plot of the function i( H) from equation (12). 
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Figure 8.— Experimental apparatus. 
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Figure 9*— Detail of short— shock— tube apparatus 
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Figure 10.- Typical voltage record obtained when shock is reflected from 

end of tube. 
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Figure 11.- Theoretical reflected-pressure ratios and measured crystal output 

for 1 -foot tube. 








Figure 13.- Comparison of theoretical and measured reflected-pressure ratios for 

11. 3 -foot tube. 
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Figure 14. - Experimental results of hydraulic tests compared with theory. 



